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I. INTRODUCTION 

Since the original suggestion by Andreev and Lif- 
shitz 1^, Chester [|| and Leggett Q that quantum crys- 
tals such as ^He might exhibit a phase where superfluidity 
coexists with crystalline order the topic of supersohds has 
received considerable theoretical attention. Early theo- 
retical work (for a review see e.g. Ref. focused on the 
possible implications of large zero-point vibrations in es- 
tablishing a supersolid phase in quantum crystals. It was 
argued by Andreev and Lifshitz 0| that these quantum 
effects might be sufficient to delocalize either impurities 
and/or zero-point vacancies. Such a system would form a 
weakly interacting Bose gas and therefore be a beautiful 
example of a system which shows Bose-Einstein conden- 
sation. Upon exploiting a mapping between hard-core 
lattice gas models and spin- 1/2 Heisenberg models it was 
shown that within a mean-field approximation such a su- 
persolid phase exists for systems with a finite range of 
interactions between the bosons [||J^. 

Unfortunately, there is not yet clear experimental evi- 
dence for such a phase. The most promising candidates 
for experimental systems, where supersolid order might 
be observed, are Josephson junction arrays and "^He films 
on substrates in two dimensions. There have also been 
experimental searches for supersolid order in bulk *He. 
In ultrasound studies 0] of highly purified solid ''He a re- 
cently observed resonance phenomenon was interpreted 
to be consistent with the presence of a supersolid induced 
by zero-point vacancies. But the experimental situation 
still remains controversial and additional evidence is 
needed to unambiguously prove the existence of such a 
phase. Artificially fabricated Josephson junction arrays 
are a particular interesting system for the observation of 
exotic phases in quantum systems. Finite-range interac- 
tions and frustrations present in these arrays give rise to 
a rather rich structure of the phase diagram. 

Partly motivated by this controversial experimental 
situation, recent work has concentrated on mapping out 



the phase diagram of strong interacting clean bosonic sys- 
tems mainly through numerical and mean-field approxi- 
mation p-p^ . The starting point of these investigations 
are lattice models of interacting bosons with the following 
Bose-Hubbard Hamiltonian, 



n 



E 



fi^^Ui + ^^niUijnj. (1.1) 



Here a^, al are boson annihilation and creation opera- 
tors at site i, and rii = a|ai. The hopping integral t sets 
the energy scale, n is the chemical potential, and Uij de- 
notes the interaction between the bosons. Bosons only 
hop to adjacent places {{i, j)). It is found that for nearest 
neighbor interaction supersolid order at half-filling is fa- 
vored at large nearest neighbor interaction Ui and small 
on-site interaction Uu = Uq, but becomes suppressed in 



the hard-core limit |14|. In the latter case next-nearest 
neighbor interactions become necessary for supersolid or- 
der to exist Besides the quantum ground states at 
half filling, which may be either a checkerboard or stripe 
domain density wave, additional Mott-insulating phases 
with a density wave commensurable with the underlying 
lattice appear. Each of these phases seems to have an as- 
sociated supersolid phase in which the particular charge- 
density order coexists with off-diagonal long range order. 
As we will show here, however, there are intermediate 
supersolid phases which have no Mott-insulating partner 
in the phase diagram. 

In the present work we report on a mean-field analy- 
sis of the hard-core next-nearest neighbor Bose-Hubbard 
model supplemented by a linear spin-wave analysis. Our 
purpose is to map out the phase diagram over the whole 
range of parameters of nearest and next-nearest interac- 
tion. This extends previous work by Bruder et al. 
which has been restricted to a small value of the next- 
nearest neighbor interaction, where the ground state pos- 
sesses checkerboard symmetry. Furthermore, we deter- 
mine the order and universality class of the phase transi- 
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tions and identify various interesting multicritical points. 

We focus on strong on-site interaction Uu = Uq, i.e. a 
hard-core approximation where the particle number per 
site takes only values between the integer values n and 
n + 1 (in the vicinity of half- filling) . Then the model can 
be transformed in a spin-1/2 XXZ model |^ 

(«J> (id) 

SfS^^^hY^St, (1.2) 

where we have dropped a constant energy offset. The 
mapping is performed by identifying a\ — Sf + iSf, 
rii = Sf-\-\, J = t and h = fJ.—J2i Uoi- Here (i, j) denotes 
the nearest neighbors and the next-nearest neigh- 

bors. After a Fourier transform, ~ '^J2q^^'^^''^q^ 
the Hamiltonian can be written as 



H — Jq [S^S^g + S^S^g) 

q - 

+ {Ulq + U2q)SlS'_^ 

with the interaction terms 

Jq = ZlJ^q, 
Ulq = -ZlUl-fq, 
U2q = -Z2U2Vq- 

The structure factors 



Sq = \l^{a^q + al), 



NhSl^^ (1.3) 



(1.4a) 
(1.4b) 
(1.4c) 



7g 



— e*^"^', and = — 

Zl Z2 



/Gn.ii. 



E 

/^n. 11.11. 



(1.5) 



depend on the lattice structure and are defined as sums 
over the nearest neighbor sites zi and next-nearest neigh- 
bor sites 02, respectively. 

In this paper we use classical ground state analysis 
and linear spin-wave theory to explore the phase diagram 
of the hard-core Bose-Hubbard system with nearest and 
next- nearest neighbor "charge" -interaction, Uij, and a 
nearest neighbor hopping term, J . We focus our study 
on the bipartite square lattice. For a discussion of frus- 
trated two-dimensional lattices, like the triangular and 
Kagome lattice, we refer the reader to a recent paper by 
Murthy et al. |l|l. 

Linear spin-wave theory is implemented in its stan- 
dard form. First, one calculates the classical ground state 
energy and the corresponding spin configuration, which 
frequently may be described by some ordering wave vec- 
tor q. Next, local rotations of the spins are performed, 
Sf Sf , such that the 2;-component points along the 
direction of the classical spin configuration. Upon per- 
forming a Holstein-Primakoff transformation, which in 
Fourier space is given by (for those ground states which 
can be described by a single wave vector q) 



Sg — 



^SSq^q 4-9-9~"P' 



(1.6a) 

(1.6b) 
(1.6c) 



a linearized spin- wave Hamiltonian Tigw: describing the 
elementary excitation around the classical ground state, 
is obtained. 

The outline of the paper is as follows. In section II we 
study the phase diagram of the hard-core Bose-Hubbard 
model with nearest and next-nearest neighbor interac- 
tion. The analysis is done by employing the mapping 
to the spin-1/2 XXZ model and using a mean-field anal- 
ysis supplemented by a linear spin-wave analysis. This 
allows us to determine the ground state spin configura- 
tions in the spin-1/2 XXZ model and analyze its stabil- 
ity and soft-modes by linear spin-wave theory. The re- 
sulting phase diagrams are given in section III together 
with a discussion of the universality classes of the criti- 
cal phenomena at the various phase boundaries. We also 
identify a series of multicritical points and discuss their 
multicritical behavior. Our results and conclusions are 
given in section IV. Finally, some technical details about 
the linear spin- wave analysis of the Mott insulating 3/4 
phase and the checkerboard and stripe domain supersolid 
phases are deferred to the appendices. 



II. MEAN-FIELD AND LINEAR SPIN- WAVE 
THEORY 



In this section we analyze the phase diagram of the 
next-nearest neighbor hard-core Bose-Hubbard model, 
using a mean-field analysis supplemented by a linear 
spin- wave analysis. However, before discussing the ef- 
fect of next-nearest neighbor interactions let us shortly 
summarize the topology of the phase diagram when the 
interaction Uij is restricted to nearest neighbors. The 
phase diagram of the spin-1/2 XXZ Heisenberg model be- 
comes particularly simple |^ (s. Fig. |l|). There is a "half- 
integer" lobe (Neel phase) centered around zero magnetic 
field {J < Ui). In the original Bose-Hubbard model this 
corresponds to the half-filling Mott-insulating phase with 
a checkerboard charge density wave. For increasing mag- 
netic field (away from half-filling) this phase becomes 
unstable to a canted ferromagnetic state, i.e. the corre- 
sponding bosonic system becomes superfluid. Finally, for 
strong fields the system goes into a paramagnetic state 
with magnetization pointing in the field direction. Such 
a state has uniform density of bosons and is insulating 
again. 
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FIG. 1. Mean-field phase diagram of a two-dimensional 
hard-core Bose-Hubbard with nearest neighbor interaction 
Ui = 1 and S — 1/2. The Neel, canted-ferromagnetic and 
paramagnetic phase in the spin- 1/2 XXZ model correspond to 
a Mott insulating phase with a checkerboard charge ordering, 
superfluid and a insulating phase with a homogeneous charge 
density, respectively. Second and first order phase transitions 
are indicated by solid and dashed lines, respectively. 



action. For that we need the wave vector dependence of 
the interaction parameters in the linear spin- wave Hamil- 
tonian. On a square lattice (xy-plane) there are four 
nearest and four next-nearest neighbors {zi — Z2 = i) 
and the structure factors are given by (a = 1) 
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(cos 



'1y) 



(2.1) 



(cos {qx + qy) + cos {qx - Qy)) = cosq^ cosqy. (2.2) 



We proceed in section [I A by a soft-mode analysis of 
the collinear phase found in the nearest neighbor model. 
This allows us to identify the stability boundaries of these 
ground states. Together with a soft-mode analysis of 
the oth er M ott-insulating phases with quarter-filling in 
section |II Bi this gives the regions in parameter space 



where intermediate phases (supersolid phases) are to be 
expected. Using an Ansatz with a general spi n co nfigu- 
ration on a 2 x 2 unit cell leads us in section |ll C| to an 
identification of the symmetry of the order parameters in 
the corresponding non-collinear phases. 



Taking into account longer-range interactions, the 
phase diagram acquires further structure. In addition 
to the phases discussed above, one finds more general in- 
sulating phases and superconducting phases in which the 
spins arrange in a 2 x 2 unit cell. In particular one can 
identify regions in parameter space where coexistence of 
long-range order in the a;?/-direction and staggered mag- 
netization in the z-direction appear. This implies that 
the corresponding bosonic system would be a supersolid, 
since there is coexistence of crystalline order and super- 
fluidity. A recent analysis of the classical ground state 
energies by Bruder et al. j|] identifies two different su- 
persolid phases, which are characterized by two or three 
angles in the 2x2 unit cell. However, their analysis has 
been restricted to small values of the next-nearest neigh- 
bor interaction, where the ground state possesses checker- 
board symmetry. Here we extend the mean-field analysis 
of Ref. to the full range of the parameters, where the 
ground state of the Bose-Hubbard model at half-filling 
possesses stripe symmetry. In addition we supplement 
the analysis of the classical ground state energy by a lin- 
ear spin-wave analysis. In performing a spin-wave anal- 
ysis we can check whether the ground state found by 
energy considerations is stable. In general, one has to 
minimize within a small symmetry class (an Ansatz), i.e. 
we make an Ansatz for the number of sublattices we need. 
From stability investigations we can confirm the ground 
state analysis, and we can evaluate the boundaries to the 
adjacent phases which occur as a soft-mode. In addi- 
tion, the wave vector of the soft-mode gives a hint of the 
adjacent phase structure. Moreover, we obtain the con- 
tribution of the quantum fluctuations to the reduction of 
the true ground state energy and the order parameter. 

Now we turn to a linear spin-wave analysis of the 
model, taking into account next-nearest neighbor inter- 



A. Collinear phases 

In this section we investigate the collinear phases, i.e. 
those phases in which the spins are parallel. These 
are the antiferromagnetic phases, namely the Neel and 
stripe phase, the paramagnetic phase and the canted- 
ferromagnetic phase. 



1. Antiferromagnetic phases: 
Neel and stripe domain ground states 

For vanishing fields, i.e. half-filling, we have a Neel, 
stripe or a ferromagnetic phase depending on the rela- 
tive strengths of the interaction parameters Ui , U2 and 
J. In the antiferromagnetic phases (Neel or stripe) the 
spins are oriented along the z-axis, which corresponds to 
an Ising-like anisotropy. Anticipating such a commensu- 
rate structure with ordering wave vector q, which is qo = 
7r(l, 1) for the Neel and qi = 7r(0, 1) or q2 — 7r(l,0) for 
the stripe phase, a Holstein-Primakoff transform yields 
the spin- wave Hamiltonian 



E 



9 L 



The coefficients of the Hamiltonian read 



Ag = S[2{Ulg + U2q) -Jq-Jc 



<2+9Ji 



B„ 



q+q 



Jq)- 



(2.3) 



(2.4a) 
(2.4b) 



For the Neel phase there are only two ground states, 
whereas four ground states with the same classical energy 
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exist for the stripe phase^ Diagonalization of the spin- with from Eq. 2.t. From the latter we conclude that 

this reduction is independent from the field, too. Calcu- 



wave Hamiltonian, Eq. 2.3, gives the classical ground 
state energies 



'ANS^{Ui - U2), Neel 
'ANS'^U2, Stripe 



(2.5) 



and the corresponding two branches of the spin-wave 



spectra — eq ± h with 



= 85 ■ 



{U2 - ^cosq^Y - j.Pcos'^qy 



(2.6) 



where we have picked the ordering wave vector qi for 
the stripe domain phase. If one takes (72 instead, then 
the corr esponding excitation spectrum is obtained from 
Eq. 2^ by interchanging qx with qy. 

Thus, when comparing the ground state energies, there 
is a transition from the Neel state to the stripe state at 



U2 = Ui/2. 



(2.7) 



Since the two ground states have different symmetry the 
transition is discontinuous (there is no intermediate spin 
orientation in a mean field approximation and from spin- 
wave analysis). The domain of stability of each of these 
classical ground states can be estimated from a soft- 
mode analysis within linear spin-wave theory. The en- 
ergy needed to excite a spin- wav e at wave vector q in the 
collinear phases is given by Eq. |2.6| . By increasing the 
field h a point can be reached where the lower branch of 
the excitation energy vanishes at a certain wave vector q. 
As a consequence, the spin- wave at this particular wave 
vector becomes soft, the classical ground state under con- 
sideration becomes unstable, and a phase transition to a 
new phase takes place. Here the minimum in the lower 
branch of the excitation energies is reached at the Bril- 
louin zone center {q — 0) for both ground states. The 
corresponding upper stability boundaries for the antifer- 
romagnetic (af ) zero-field ground states are given by 
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_ / 9.S^{Ui~U2Y - J^, Neel 



SS^U2{U2 - J), 



Stripe 



(2. 



For fields h lower than h'^{^^ 



the Neel and stripe domain 
phases are stable against quantum fluctuations. Note 
that in the stripe phase the parameter Ui does not occur 
explicitly in the ground state energy and in the spin- wave 
spectrum. 

For h < h™^''^ the dispersion relation shows an energy 
gap for both phases because there is no continuous de- 
generacy of the ground state. Thus, these phases have 
long-range order even at finite temperatures. At T = 
the order parameter (staggered magnetization) is reduced 
by quantum fluctuations 



-y 



2NS 



(2.9) 



lating the magnetization we readily see that there is no 
contribution from quantum fluctuations; thus, the mag- 
netization vanishes in the whole lobe. This corresponds 
to a constant mean particle number (n) ~ {S^}, which is 
a signature of a Mott insulator. The tip of the lobes of 
the antiferromagnetic phases can be found from the zeros 
of the critical fields, Eq. 2.8. The corresponding limits 



for the hopping integrals are 



J < 



Ui — U2, for checkerboard symmetry 
U2, for stripe symmetry. 



(2.10) 



2. Paramagnetic phase 

In the paramagnetic phase all spins are oriented along 
the field direction (z-direction). Applying the Holstein- 
Primakoff transform we get the following spin-wave 
Hamiltonian 



Hsw — Eg + ^ <^qa\aq , 



(2.11) 



with the classical ground state energy Eg and the spin- 
wave spectrum 



Eg = ANS'^iUi + U2) - NSh, 
eg^h-2S{Uw + U2o + Jq). 



(2.12) 
(2.13) 



Stability of the ground state requires a positive excitation 
spectrum which gives the lower bound of the paramag- 
netic phase (soft-mode for q = 0) 

= hc = SS{Ui + U2 + J). (2.14) 



3. Canted-ferromagnetic phase 

For vanishing field and large values of J the spins align 
ferromagnetically in the plane. Due to the rotation sym- 
metry in the xy-plane there is a Goldstone mode. For 
an infinitesimal field perpendicular to the plane this fer- 
romagnetic phase becomes unstable and changes to a 
canted phase where the spins orient ferromagnetically to- 
wards the field direction. By minimizing the energy we 
get a relation between the canting angle and the magnetic 
field 



h — hr sint 



(2.15) 



Here 9 denotes the angle between the plane and the spin 
direction. The resulting ground state energy can be writ- 
ten as 
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Eg = -ANS^J - NS— , 



(2.16) 



and the dispersion relation for the canted-ferromagnetic 
state is 



Ja 



Jq ~ Jq 



(2.17) 



By comparing the classical ground state energy of the 
paramagnetic and the canted-ferromagnetic phase one 
finds that the transition to the paramagnetic phase takes 
place at h = hg. This agrees with the lower stability 



boundary, Eq. 2.14, obtained from a soft-mode analysis 
of the paramagnetic phase. At the transition the cant- 
ing angle continuously goes to tt/2. Thus, in a mean- 
field approximation the corresponding phase transition 
a.t h = he is of second order and belongs to the D — 3 
classical XY universality class. 

The lower field boundary can again be obtained from 
a soft-mode analysis of the excitation spectrum. One 
finds that depending on the relative magnitude of the 
next-nearest and nearest neighbor interaction there is a 
soft-mode either at the ordering wave- vector of the Neel 
state q = qo 01 of the stripe domain state q = qi, q2- 
This gives 



Ncol 



Ul~U2-J 



C/1-C/2 + J' 

%S{Ui + U2) ~ 16S' 



for go, 
V2 



C/1-C/2 



J, 



and 



U2~J 



U2 



for qi, 92- 



U2 



(2.18a) 
(2.18b) 

(2.19a) 
(2.19b) 



At these boundaries the canted-ferromagnetic phase be- 
comes unstable to an intermediate phase with a checker- 
board or stripe-li ke sp in configuration which will be dis- 
cussed in section II C . Here we added the expansion for 
small J . 



B. The Mott-insulating 3/4 and 1/4 lobe 

In the preceding section we have seen that there must 
be intermediate phases with a more complicated spin con- 
figuration. As shown in previous work there exist 
phases with non-integer filling, namely 3/4 and 1/4, i.e. 3 
up and 1 down spin and vice versa for a reversed magnetic 
field. The spins have no component in the plane. Such 
a ground state on a square lattice corresponds to a Jour- 
sublattice system for which the ground state is uniquely 
defined. By symmetry arguments the down spin (3/4 



lobe) can be at each position without changing the en- 
ergy. The classical ground state of this configuration is 
given by 



En = -IrNSh, 



(2.20) 



which is independent of the exchange energy J, Ui and 
C/2. Performing a spin- wave analysis we can calculate the 
excitation energy of four branches. From stability condi- 
tions, i.e. positiveness of the excitation energies, we can 
deduce various equations. For J — the Hamiltonian is 
already diagonal and the four branches of the dispersion 
relation are given by (s. appendix) 



9 
^(2) 

g 

.(3) 



h 
h 



= h 



8SU2, 
8S{U2 
h + 8S{U2 
- 8SU2. 



(2.21a) 
(2.21b) 
(2.21c) 
(2.21d) 



These excitations are independent of the wave vector 
q. From these equations we get an upper and two 
lower bounds of the 3/4- lobe, depending on the relative 
strength of Ui and U2: 



8SiU. + U2)>h>[ll\^^^^ 



U2 



Ui > 2U2 
Ui < 2U2 



(2.22) 



The 3/4 phase vanishes for vanishing 1/2- The location of 
the lower bounds are different, depending on whether the 
zero-field ground state possesses checkerboard or stripe 
domain symmetry. For finite J one has to solve the full 
dynamical matrix. Vanishing of this spectrum at g = 
leads to the following condition (see appendix ^) 



^ - 8S{2Ui - 
i8SfU2{U^ 



U2)h^ 



{8Sf{Ul 
J') = 0, 



2UiU2)h 

(2.23) 



which determines three solutions. The fourth solution 
defines the boundary h = 8SU2, which corresponds to 
e*-"*) and is independent of J. To see which solutions are 
relevant we expand them for small J values: 



8SU2 



8SU2 



Ui{2U2~Ui) 



J' 



~85([/2-C/i) + 
~ 8S{JJi + f/2) - 



AS 



Ui - 2U2 



-J' 



45 

IT, 



(2.24a) 

(2.24b) 

(2.24c) 
(2.24d) 



From these expansions we see that /i'^' defines the upper 
boundary and /i^^^ and ft.^^-* the lower boundaries for the 
stripe and the Neel ground states respectively, /i^^^ is ir- 
relevant. As will be shown in the next section, the lower 
boundaries h'^^^ and /i'^' define a continuous transition 
to a non-coUinear three-sublattice structure (SS2). 
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In comparison with the phase boundaries for the 



colhnear spin phases (Eq. 2.10) we readily see that the 
antiferromagnetic phases only meet the 3/4 phase at 
J = 0. However, the upper boundary needs some more 
investigation: the phase boundaries from the canted- 
ferromagnetic to a non-coUinear phase (Eqs. 2.18b and 
2.19t) show that the 3/4 phase intrudes (for U2 < Ui) 
in the canted ferromagnetic phase. The solution of this 
seemingly paradoxical situation is that there is a discon- 
tinuous transition from 3/4 to the canted-ferromagnetic 
phase, (for low J) which is given by the equality of their 
ground state energies: 

h,3 = 45([/i + U2 + J) 



+4.S^{Ui + U2 + J){Ui + U2- 3 J) (2.25) 



8(C/i 



C/2) - 



85* 



U1 + U2 



Thus, for U2 < Ui this discontinuous transition is pref- 
ered and the upper boundary (/i'^'^^) defines a continuous 
transition to a supersohd phase (SS2) only for U2 > Ui. 



C. Non-collinear phases: supersolids 

The stability analysis in the preceding sections has 
shown that there must be several non-collinear phases 
between the paramagnetic and the Neel/stripe phase in 
which the spins have both a ferromagnetic and an anti- 
ferromagnetic orientation. In order to describe these in- 
termediate phases we use an Ansatz with a general spin 
configuration with four different angles ai between the 
z-axis and the xy-plane, i = l,---,4, ina2x2 unit cell 
(see Fig. I). 




FIG. 2. General spin configuration with four different an- 
gles Qfi, i = 1, • ■ ■ , 4, in a 2 X 2 unit cell. 

Due to the planar symmetry only these angles are in- 
dependent. In such a general four-sublattice model the 
ground state energy is given by 



cos ai cos 



{i,j)en.N. 



-2U2 E 



Uj — Ui sin ai sin aj 

{i,j)en.N. 
4 



sma, sma., 



1 



hS^sina^, (2.26) 



where the sums run over the nearest neighbor (X^n n ) 
and next-nearest neighbor {J2n n n ^ sites within the 2x2 
unit cell. 



1. Two-sublattice supersohd phases: SSI and SSI* phase 



Starting from a general spin configuration with four 
different angles on a 2 x 2 unit cell there are two different 
possibilities for a general two-sublattice structure. Cor- 
responding to the Neel and the stripe domain zero-field 
structure one may have spin configurations with checker- 
board («! = Q!3 = a, a2 = Q!4 = (3) and stripe sym- 
metry (cki — a2 — a, as = a4 = P or ai = a4 — a, 
a2 — as — (3), which we term SSI and SSI*, respectively 
(see Fig. I). 



FIG. 3. Sketch of the checkerboard and stripe symmetry. 
The left figure corresponds to SSI and the two right ones to 
SSI*. The four squares represent the four spin angles 01-04, 
where the same color means the same angle. 

Using such a two-sublattice Ansatz the corresponding 
ground state spin configuration is obtained by minimiz- 
ing the classical energy, Eq. 



2.26. This gives the two 



angles a and /? as a function of the magnetic field h and 
the interaction parameters. By comparing the resulting 
ground state energy with the corresponding ground state 
energies of the Neel, canted-ferromagnetic and 3/4 phase 
upper and lower bounds for the stability of the super- 
solid phase can be obtained. In order to show that these 
bounds are actually phase boundaries, they must be com- 
pared with the results obtained from a soft-mode analy- 
sis. Furthermore, it must be determined whether there 
are additional intermediate phases with a spin structure 
more general than the two-sublattice Ansatz discussed in 
this subsection. 

After those general remarks, let us now discuss the 
explicit results from such a mean-field analysis in the 
various parameter regions. 

(i) For Ui > 2U2 the checkerboard structure is ener- 
getically preferred, i.e. we expect a phase transition from 
the Neel phase to the intermediate phase SSI. With the 
introduction of sum and difference angles 



7 
<5 



Q + /3 

a~ (3 



(2.27) 
(2.28) 



minimization of the classical energy leads to a spin con- 
figuration described by the angles 
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7 = TTT^ttta/ ^' ^rr -Ah Eo), (2.29) 



I6SU2 \ Ul^U2 + J 



2r 1 U1-U2 + J 

The ground state energy reads 
N 



{h - Eo). (2.30) 



En 



where 



32U2 



[{h - Eof + 2{SSfU2{Ui - U2)\ , (2.31) 



So = BS^{Ui-U2- J){Ui^U2 + J). 



(2.32) 



Comparing wi th the ground state energy of the Neel 
phase, Eq. 2^, one gets the lower bound of stabihty of 
the SSI phase 



= So = 85V(C/i -U2- J){Ui -U2 + J), (2.33) 

which is identical to the upper stability boundary ob- 
tained from a soft-mode analysis of the Neel phase, 
Eq. 2.8. Therefore, one may conclude that this is the 



actual phase boundary between the Neel and SSI phase, 
and the corresponding phase transition is continuous (at 
the mean- field level). The high- field boundary to the 
canted-ferromagnetic phase is obtained by comparing 
with the corresponding classical gro und s tate energy of 
the canted-ferromagnetic phase, Eq. 2.16| , or more con- 
veniently by the condition a = (3. One finds 



hc2 = K 



Ul~U2~J 

U1-U2 + J' 



(2.34) 



which agrees with the lower stability bound obtained 
from a sof t-mode analysis of the canted-ferromagnetic 
phase, Eq. 2.18t| . Hence, this is again an actual phase 
boundary marking a continuous transition between the 
canted-ferromagnetic and the supersolid SSI phase. Of 
course, this conclusion is valid only if the transition is not 
pre-emptied by a transition of the canted-ferromagnetic 
phase to another phase which is lower in energy than 
the SSI phase. As already seen in the previous sec- 
tion there is a discontinuous transition from 3/4 phase 
to the canted-ferromagnetic phase for low J values. The 
two phase-boundaries meet at the multicritical point mci 
(s. Figs. 1^ and ^) determined by the implicit equation 
hc2 = hc3. Finally, by comparing the classical ground 
state energies of the SSI and the 3/4 phase one obtains 
the stability boundaries 



So + 8SU2 

±^16SU2{En + iS{3U2 - 2Ui)). 



(2.35) 



Due to different symmetries (two-sublattice for the SSI 
and three-sublattice for the 3/4 phase) this transition 
is discontinuous. This transition can appear only for J 
larger than the values at the multicritical point mci. It 
meets the phase boundaries hc2 and /ic3 at mci and hd 



at J — for a certain region of the parameters Ui and 
U2- Without a stability analysis of the SSI phase, it is at 
this point not possible to conclude that the line defined 
by Eq. is an actual phase boundary. As we will see 
in the following sections the conclusion depends on the 
actual choice of the parameters Ui, U2 and J. 

(ii) For Ui < 2U2 the low-field phase is the antiferro- 
magnetic stripe phase. The analog intermediate phase 
(SSI*) has a ground state energy of 



F ^ 
' 16(C/i + J) 



[ih-E*f + i8SfU2iU, + J)] 

(2.36) 



with 



E*o = 8S^U2{U2 



J). 



(2.37) 



Here the lower boundary, h = E*, equals Eq. 2.8 and the 
upper boundary is 



K2 = 



U2~J 



(2.38) 



which agrees with Eq. 2.19h| . In the limit of small J 
values this transition is given by 



8S{Ui + U2 



■AS 



U2-U. S{Ui + 5U2) 



U2 



-J 



J' 



(2.39) 



Thus, the slope is positive for U2 > Ui, in contrast to 
all the other transition lines at this point (except h^^^). 
We will see later that this influences the topology of the 
phase diagram. 

The relation between the angles and the field is in anal- 
ogy to the checkerboard case 



sin^ 7 



cos^ S 



1 



U2-J 



8S{Ui + J) V U2 

1 r~u7 



{h~E*), 



8S{Ui + J)\IU2~J 



;ih~E*,). 



(2.40) 



(2.41) 



Performing a stability analysis of both phases via linear 
spin-wave theory (see appendix ^ shows the occurrence 
of another phase, a three-sublattice phase SS2. The ex- 
citation spectrum for the SSI phase has a soft-mode at 
q = 7r(l,0) and q = 7r(0, 1), whereas the SSI* has a 
soft-mode at q — 7r(l,l). The transition line has been 
determined numerically from the expressions in the ap- 
pendix. The soft-mode indicates a continuous transition. 
However, the analysis of the SS2 phase in the next section 
shows that this is only partially true. 
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2. Three-sublattice supersolid phase: SS2 and SS2* 

From spin-wave considerations in the last section we 
encountered another phase which has to be more comph- 
cated than the general two-sublattice structure. There- 
fore we are led to consider a three-sublattice structure 
(SS2 phase) and calculate the phase boundaries to the 
other phases discussed above. We take the following pa- 
rameterization of the angles (s. Fig. ^ 

as = ai = a, "2 = + Ci "4 = ?7 - C 7 (2.42) 
which has partially a checkerboard symmetry. 



FIG. 4. Sketch of the checkerboard and stripe symmetry. 
The left figure corresponds to SS2 and the right one to SS2* 
which is not stable in a four-sublattice model. The four 
squares represent the four spin angles 01-04, where the same 
color means the same angle. 

In sertin g in the general four-sublattice expression 
(Eq. 2.26 ) for the ground state we can evaluate the mean- 
field equations for the angles and the resulting ground 
state energy. Assuming a continuous transition from the 
SS2 phase in the 3/4 phase we recover the result al ready 
found for the 3/4 phase by spin- wave analysis, Eq. 2.23 . 
This equation, as already discussed, defines only partially 
a continuous transition. Searching for a continuous tran- 
sition to the SSI phase by comparing the corresponding 
ground state energies we obtain the following equation 



hU2Eo sin (7 - 5) = 8S[Ui{Ui - U2) - J^]h 
+8S{U^ ~ Ul -f J2) , 



(2.43) 



which defines the same boundary as found numerically 
by the spin-wave analysis of the SSI phase. The an- 
gles 7 and 5 are given by the solution for the SSI phase 
(Eq. |2.30 ). Investigation of the ground state by numer- 



ical minimization (simulated annealing , see next subsec- 
tion) of the ground state energy, Eq. 2.26 , shows that 
there are regions where the transition is continuous and 
regions where the transition is discontinuous. We also 
find that there is no continuous transition to the SSI* 
phase possible due to the different symmetry. 

A corresponding SS2* spin structure with stripe sym- 
metry can be parameterized by (s. Fig. ^ 

a2 = ai = a, — rj-{- ai — rj — Q. (2.44) 

Numerical investigation shows that this phase seems to 
be absent in the whole parameter range. Within our nu- 
merical resolution we could not find any significance for 
it. From symmetry considerations we expect that there 



is no phase of that kind. There are only three different 
two-sublattice wave- vectors 90,(71 and (72 with which we 
cannot construct the symmetry of a SS2* phase. 



3. Four-sublattice supersolid phase: SS3 

In order to study the general four-sublattice structure 
we have employed two different kinds of numerical meth- 
ods. Using a simulated annealing code we have deter- 
mined the ground state spin configurations on the 2x2 
unit cell; these simulations were crosschecked by a stan- 
dard library for solving nonlinear equations. We find 
that the general four-sublattice structure appears only for 
Ui < 2U2, i.e. for the stripe phase ground state. Thus 
we find an intermediate supersolid phase (SS3), which 
interpolates between the SSI* and SS2 phase such that 
both transitions become continuous. The resulting phase 
diagrams are shown in the next section for various values 
of the interaction parameters Ui, U2 and J. 



III. PHASE DIAGRAMS, MULTICRITICAL 
POINTS AND UNIVERSALITY CLASSES 

In this section we summarize our results for the phase 
diagram of the hard-core Bose-Hubbard model with next- 
nearest neighbor interaction, as obtained from the lin- 
ear spin-wave analysis. Furthermore, we identify various 
multicritical points and comment on the possible univer- 
sality classes of the phase transitions at the second order 
phase boundaries. 

A number of different phase diagram topologies are 
possible depending on the relative magnitude of the near- 
est and next-nearest neighbor interaction. Three differ- 
ent classes may be distinguished. For small values of 
the next-nearest neighbor interaction, U2 < Ui/2, the 
zero-field ground state is characterized by checkerboard 
symmetry. In an intermediate range, Ui/2 < U2 < Ui, 
the zero-field ground state switches to stripe symme- 
try and there is strong competition between nearest and 
next-nearest neighbor interaction. Finally, the third class 
of phase diagram topologies is obtained when the next- 
nearest neighbor interaction is dominant, U2 > Ui. 



A. Type I topology: Checkerboard ground state and 
small next- nearest neighbor interaction (U2 < Ui/2) 

For relatively small next-nearest neighbor interactions, 
U2 < Ui/2, the classical ground state at zero magnetic 
field is the Neel state. If the interaction is restricted to 
nearest neighbors {U2 = 0) there are only three phases: 
the Neel, the canted-ferromagnetic and the paramagnetic 
phase (s. Fig. |l|) . Taking into account longer-range inter- 
actions, i.e. for finite U2, an additional insulating phase 
with three quarter filhng (3/4 phase) and two supersohd 
phases, SSI and SS2, appear 
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FIG. 5. Part of the phase diagram for the hard-core 
Bose-Hubbard model for Ux = \ and U2 = 0.1 (cf. |||). Sohd 
hnes indicate continuous phase transitions (to mean-field 
level), whereas first order phase transitions are given by 
dashed lines. At finite J there are four multicritical (mc) 
points: (i) two critical end points mci and mc2, (ii) a tricriti- 
cal point on the tip of the supersolid SS2 phase (marked by a 
star (*)), and (iii) a multicritical point mcs on the tip of the 
Neel phase, which is not shown here. 
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FIG. 6. The same as in Fig. ^ for interaction parameters 
C/i = 1 and 172 = 0.2. 
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FIG. 7. The same as in Fig. ^ for interaction parameters 
= 1 and (72 = 0.45. 



In Figs. |^-|^ we have shown the phase diagrams re- 
sulting from ground-state calculations and linear spin- 
wave analysis for U\ — \ and Ui — 0.1, C/2 — 0.2, and 
C/2 = 0.45, respectively. The characteristic features of 
the mean- field phase diagram are as follows. 

(I) Multicritical points: At finite values of the hopping 
matrix element there are four multicritical (mc) points 
mci, mc2, mc3 and the tricritical point on the tip of the 
SS2 phase (marked by a star in Figs. Since for a 

finite value of the next-nearest neighbor interaction U2 
there is an intermediate supersolid phase intervening the 
Neel and the canted-ferromagnetic phase, there is a mul- 
ticritical point mc3 at the tip of the Neel phase where 
two second order phase boundaries meet each other with 
a common tangent. For U2/U1 < 0.183 the second order 
line separating the canted-ferromagnetic (i.e. superfluid) 
from the supersolid SSI phase intersects and is truncated 
by a first order phase boundary to the commensurate 3/4 
solid phase. The corresponding multicritical point mci 
is a quantum critical end point. The "spectator phase" 
is the uncritical Mott insulating 3/4 phase. The second 
quantum critical end point mC2 differs from mci in one 
important aspect. This time the "spectator phase" is 
the supersolid SSI phase characterized by the simultane- 
ous presence of crystalline order and superfluidity. Hence 
the "spectator phase" is critical since it shares the super- 
fluid order and the corresponding Goldstone mode for the 
phase of the superfluid order parameter with one of the 
phases (the SS2 phase) associated with the critical line. 
Fisher et al. ||l^,|l^ have shown that the singular behav- 
ior at the critical end point also engenders new singular- 
ities in the flrst order phase boundary itself. Therefore 
one may speculate that the presence of the critical end 
point may also affect the "spectating" supersolid phase 
and induce "non Bose-liquid" behavior. 

Upon increasing the relative magnitude of the next- 
nearest neighbor interaction the two critical end points 
mci and mc2 approach each other and for U2/U1 w 0.183 
they merge into a higher order multicritical point. The 
critical values at this point are Jc — 0.365 and hcnt — 
3.827. From a higher dimensional perspective, where the 
ratio U2/U1 is added as an additional axis in the phase 
diagram, this point corresponds to the crossing point of 
two multicritical lines. In the projection of this higher 
dimensional phase diagram, mci and mc2 switch their 
relative position for U2/U1 > 0.183: now at mc2 the 
3/4, the canted-ferromagnetic and the SS2, and at mci 
canted-ferromagnetic, SSI and SS2 meet. 

There is also a tricritical point on the tip of the SS2 
phase where the first order line separating the two super- 
solid phases terminates and becomes second order. The 
presence of the tricritical point is to be expected, since 
critical end points in classical critical phenomena are in- 
timately associated with the vicinity of tricritical points 
|T7| . With increasing the next-nearest neighbor interac- 
tion the critical end point mci and the tricritical point 
approach each other. At U2 ~ Ui/2 they merge into a 
higher order multicritical point. 
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(II) Supersolid phases: There are two different super- 
solid phases, SSI and SS2, and each of them has a Mott- 
insulating partner phase with the same symmetry of the 
charge ordering, the Neel and the 3/4 phase, respectively. 
The phase transition between the two supersolid phases 
is first-order at high fields and becomes second-order at 
low fields. The tricritical point, where those two tran- 
sition lines meet is located at the tip of the SS2 phase. 
With increasing U2 the SS2 phase fills more and more 
of the parameter region occupied by the SSI phase (see 
Fig. ^1^). At Ui — 2U2 the phase diagram changes dras- 
tically because the stripe phase replaces the Neel phase. 



0.40 



0.20 



0.00 




0.05 



FIG. 8. Location of the tip of the supersolid SS2 phase 
as a function of the ratio U2/U1. The maximum of the tip 
is reached before the transition to the stripe ground state, 
namely for U2/U1 ~ 0.43 with a maximum value of J ~ 0.507. 

In Fig. [s] a plot is shown for the tip (J- value, marked by 
a star in the phase diagrams) of the SS2 lobe as a function 
of U2/U1. From the soft-mode analysis (see appendix ^ 
of the SSI phase we can derive a cubic equation which de- 
fines the tip by the collapse of two solutions. Surprisingly, 
the maximum of the tip is reached before the transition 
to the stripe ground state, namely for U2/U1 ~ 0.43 with 
a maximum value of J ~ 0.507. As it turns out the SSI* 
phase intrudes in the area of the SSI phase for U2 slightly 
smaller than 0.5 (s. Fig |l8| ). 

(Ill) Universality classes: Whereas the mean-field 
analysis in section |l^ is certainly sufficient in describing 
the overall topology of the phase diagram and the or- 
dered phases away from the transition, it breaks down 
near the phase transition where the correlation length 
diverges. The critical behavior at the phase transition 
from the SSI phase (checkerboard supersolid — XSS) to 
the canted-ferromagnetic phase (i.e. superfluid phase in 
the corresponding Bose system) has recently been stud- 
ied by renormalization group theory . It is found that 
the phase transition from the checkerboard supersolid to 
the superfiuid phase can be described by an effective ac- 
tion, where an Ising field is coupled to the superfluid 
phase fluctuations {XSS-SF universality class). It ex- 
hibits nontrivial critical behavior and appears within an 
e-expansion to be driven first order by fluctuations. How- 
ever, within a calculation directly in d = 2 dimensions 
one finds a fixed-point with "non Bose-liquid" behav- 
ior. At the transition between the two supersolid phases 



(at low fields) there are soft-modes at qi — 7r(l,0) and 
12 — ''■(0,1) (see appendix]^). Hence the boson density 
in the SS2 phase may be characterized by 

n{x) ^ uaRe (l + (Pe"""" + tPie"""" + i!2e"^^'') , (3.1) 

where no is a background boson density and the fields 4> 
and i/j represent the checkerboard and stripe domain or- 
der parameters, respectively. From the mean-field anal- 
ysis in the preceding section we know that in the SS2 
phase, the system orders along the [ll]-direction in the 
("011 V'2)-plane. This implies a particular form of the cu- 
bic anisotropy in the effective action in the SS2 phase 
{v' > 0, see below). The XY field is linearly coupled 
to the Ising field (f> describing the checkerboard order of 
the SSI phase. The corresponding action allowed by sym- 
metry reads M 



1 



(V(? 



t 



4! 



1 

+ 2 



V' 

'4! 



(3.2) 



with the control parameters t and t' measuring the dis- 
tances from the critical lines. Within the linear spin- wave 
theory in section || these critical lines coincide with those 
curves in parameter space, where the spin-wave spectrum 
has a soft-mode at the ordering wave vector of the cor- 
responding field. The parameters of the effective action 
are related to the microscopic parameters of the original 
Bose-Hubbard model, which may be worked out using a 
mean-field decoupling procedures (see e.g. the appendix 
of Ref. ) . Non-local interactions arise due to interac- 
tions with long wavelength fluctuations of the superfluid 
phase, 9, 



Si = 



(V£ 



(3.3) 



where we have used the same notation as in Ref. |1§|| . 
This part of the effective action originates from the x and 
y spin components in the spin-1/2 XXZ model, Eq. |l.2| . 
Sufficiently far from the phase boundary where the su- 
perfluid order parameter vanishes, the most relevant cou- 
pling to the spatial order parameters allowed by the time- 
reversal and U{1) symmetries is |18[ 



S2 = / dxdrl i(jdre\(t)\^ + icr'a^6l|i/>p 



(3.4) 



For the phase transition from the SSI to the SS2 phase, 
the spin-wave spectrum softens at the stripe domain or- 
dering wave vectors qi and q2, i.e. t' becomes zero. Since 
this happens in the presence of an already ordered Ising 
field (f) there are now two different cubic anisotropics. 
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the quartic term and the quadratic term, g'ipitp2, where 
g — w{(f>). Whereas the quartic cubic anisotropy is an 
irrelevant symmetry breaking field for the XY transi- 
tion the quadratic symmetry breaking field g leads 
to an uniaxial anisotropy in the [ll]-direction and hence 
to a reduction in the number of soft spin components 
from n = 2 to n = 1. This in turn implies that the 
coupling to the phase fluctuations becomes relevant to 
the asymptotic critical behavior at the SSI to SS2 phase 
transition, and the asymptotic critical behavior is in the 
XSS-SF universality class. When the checkerboard or- 
der parameter (cj)) is decreased, the quadratic symmetry 
breaking field becomes small, and there can be quite in- 
teresting crossover phenomena. One should note that 
here both the cubic and the quartic cubic anisotropy fa- 
vor spin orientation in the [ll]-direction, since v' > 0. 
Later we will encounter a case where the cubic anisotropy 
becomes negative and as a consequence of the compe- 
tition between these two anisotropics new intermediate 
supersolid phases appear. 

Interesting new quantum critical behavior is also found 
in the vicinity of the tricritical point at the tip of the 
supersolid SS2 phase. Since the tricritical point of the 
classical D — 3 n-vector model is described by classical 
exponents with logarithmic corrections, the specific heat 
exponent a becomes positive, a — 2 — Dv — 0.5. As a 
consequence the coupling to superfluid phase mode be- 
comes even more relevant on approaching the tricritical 
point. The effective action at the tricritical point reduces 
to 



5o — 



(ix(ir< — 
I 2 



r-' 



6! 



3^1 -02 



(3.5) 



Due to the positive value of a, a renormalization group 
analysis of this quantum tricritical model will yield a new 
universality class different from the D = 2> tricritical 
Ising model. We leave the corresponding analysis of the 
asymptotic critical and crossover behavior for a future 
investigation pO| |. 

Finally, the topology of the phase diagram allows for 
phase transitions between the supersolid and commensu- 
rate solid phases. At both transition the order parameter 
characterizing superfiuid order becomes zero. In passing 
from the SSI to the Neel phase the order parameter de- 
scribing checkerboard ordering changes smoothly, while a 
two-component field 'J describing superfluidity becomes 
critical. Since generically there is no particle-hole sym- 
metry (charge conjugation) symmetry along this phase 
boundary [ pT| , the effective action can contain a term 
^*dr^ which is more relevant than jcJ^^'p. Hence, the 
transition in 2 -I- 1 dimensions becomes mean-field like 
with logarithmic corrections. Similar conclusions hold 
for the transition from the SS2 to the 3/4 phase. 



B. Type II topology: Intermediate next-nearest 
neighbor interaction (Ui/2 < U2 < Ui) 



In this parameter region the zero-field ground state is 
an antiferromagnet with a striped structure. As a con- 
sequence of the strong competition between nearest and 
next-nearest neighbor interaction the phase diagram ex- 
hibits quite a rich structure with first and second order 
phase boundaries and various multicritical points. The 
characteristic features are as follows: 
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FIG. 9. Phase diagram for the hard-core Bose-Hubbard 
model with interaction parameters (7i = 1 and U2 ~ 0.55. 
The first-order and second-order lines in the phase diagram 
are depicted by dashed and solid lines, respectively. There are 
five multicritical points: (i) one critical end point, mc2, (ii) 
two bicritical points, BC and mci, (iii) one tricritical point 
along the phase boundary between the supersolid SS2 phase 
and the canted-ferromagnetic phase in the intermediate vicin- 
ity of mci (marked by a star (*)), and (iv) the multicritical 
point on the tip of the checkerboard solid phase. The small 
sliver is the SS3 phase. 
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FIG. 10. Phase diagram for the hard-core Bose-Hubbard 
model with interaction parameters Ui — 1 and U2 = 0.75. 
The bicritical point on the tip of the SS3 phase has merged 
with the bicritical mci to from a tetracritical point (QC). Oth- 
erwise the phase diagram has the same topology as Fig. ^ 
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FIG. 11. Phase diagram for the hard-core Bose-Hubbard 
model with interaction parameters Ui = 1 and U2 = 0.9. The 
phase diagram has the same topology as Fig. ^[ The tricrit- 
ical point on the phase boundary between the SS2 and the 
canted-ferromagnetic phase has shifted towards the multicrit- 
ical point mc2. 



(I) Supersolid phases: There are three supersolid 
phases. In addition to the SSI and SS2 phase we find 
an intermediate supersohd phase SS3 characterized by 
a four-sublattice structure with all four angles on the 
2x2 unit cell being different. Increasing the next- 
nearest neighbor interaction U2 from the lower bound 
U2/U1 = 0.5 a sliver of the SS3 phase appears at small 
values of the hopping matrix element J (see Fig. ||). Both 
of the transitions of the other two supersolid phases to 
this intermediate supersolid phase are continuous. This 
can also be inferred from Fig. where the variation of 
the four angles of the spin configuration are shown as a 
function of the magnetic field h for fixed hopping matrix 
element J = 0.40. 
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FIG. 12. Variation of the angles on the 2x2 unit cell char- 
acterizing the spin configuration at fixed interaction param- 
eters J = 0.40, Ui = 1 and U2 = 0.75 as a function of the 
magnetic field h. All transitions are seen to be continuous on 
the mean-field level. The boundaries of the SS3 supersolid 
phase are marked by dashed lines. 



For larger values of J the transition is still directly from 
the SS2 to the SSI* phase with a jump in the angles of 
the spin configuration, i.e. the transition is first-order. 
As the ratio of U2/U1 becomes larger the SS3 phase oc- 
cupies a larger portion in parameter space and finally (at 
U2/U1 ~ 0.75) meets with its tip the phase boundary to 
the canted-ferromagnetic phase. See the topology of the 
phase diagram in Fig. |l^, where U2/U1 = 0.75. For larger 
values of U2 the SS3 phase keeps its topology, i.e. there 
is an intermediate SS3 phase completely separating the 

552 from the SSI* phase by a series of continuous phase 
transitions. 

(II) Multicritical points: Due to the competition be- 
tween nearest neighbor and next-nearest neighbor in- 
teraction there are also several interesting multicritical 
points. Here we discuss those which appear at finite val- 
ues of J. Right at U2 = Ui/2 the critical end point mci 
and the tricritical point at the tip of the SS2 phase (com- 
pare Fig. 1^) have merged into a tricritical end point. A 
slight increase of U2 above the critical value Ui/2 implies 
that part of the phase boundary between the SS2 and 
the canted-ferromagnetic phase becomes a second order 
line. Then, mci becomes a bicritical point and a new tri- 
critical point close to mci appears (compare Fig. ^ . The 
bicritical point separates two second-order lines from the 
superfluid to two supersolid phases with distinct order- 
ing wave vectors. The SSI* phase is characterized by one 
of the stripe-domain ordering wave vectors qi = 7r(0, 1) 
or q2 = 7r(l,0); the charge ordering in the SS2 phase is 
more complicated and in addition to qi and 52 requires 
a third (checkerboard) wave vector qo — 7r(l, 1). In ad- 
dition there is a bicritical point (BC) on the tip of the 

553 phase. Upon increasing the next-nearest neighbor 
interaction the tip of the SS3 phase and hence the bicrit- 
ical point moves towards the bicritical point mci until 
both join at the critical value U2 ~ 0.75 to a tetracritical 
point (QC) |2^. At the same time the phase boundary 
connecting mci and mc2 becomes more and more of sec- 
ond order. As a consequence the tricritical point shifts 
toward the critical end point mc2 (compare Fig. [Tl| ) . Fi- 
nally, at U2 — Ui the phase boundary between the SS2 
and canted-ferromagnetic phase is completely of second 
order and the tricritical point merges with the critical 
end point to become a tricritical end point at 

J, = 5-1, h.^Ag, (3.6) 

where g = {1 + \/5)/2 denotes the golden mean. The tip 
of the 3/4 phase is located at Jtip = (2/\/27)^/^ > Jc- 
This tip appears for [/2 > 0.96. At U2 = Ui the phase 
line for smaller J values is still a discontinuous transition 
between the canted-ferromagnetic and the 3/4 phase but 
becomes a complex transition for slightly larger values (s. 
subsection C). 

(III) Universality classes: In each of the three super- 
solid phases the boson density may again be character- 
ized by 

n(x) = noRe (l + ^e**"^ + ^pie"^'"" + ^'26*'?'"=) . (3.7) 
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For U2/U1 < 0.75 and large enough J there is a first 
order hne separating the SSI* from the SS2 phase. In 
both phases the stripe domain order parameter ^p is fi- 
nite, but diflters in its ahgnment in the (i/;!, ■!/)2)-plane. 
Whereas the vector -0 points along one of the cube edges 
in the SSI* phase, it is oriented along one of the cube 
diagonals in the SS2 phase. From the mean-field anal- 
ysis in section ^ we also know that the checkerboard 
order parameter (f) is zero in the SSI* phase but finite in 
the SS2 phase. Below a certain critical value of J which 
depends on the ratio U2/U1 there is an intervening super- 
solid phase (SS3) characterized by an order parameter ip 
which rotates from the cube edges to the cube diagonals 
as one passes from the SSI* to the SS2 phase. 

The topology of this part of the phase diagram can be 
understood in te rms of a mean-field analysis of the effec- 
tive action, Eq. |3.2| . The phase boundary between the 
SSI* phase and the SS3 phase corresponds to t = and 
the phase boundary between the superfluid phase and 
both of the supersolid phases corresponds to t' = 0. The 
first order transition between the SSI* and the SS2 phase 
is a consequence of the competition between the quartic 
{v' < 0) and quadratic (w) cubic symmetry breaking 
fields in the effective action. A mean-field analysis of the 
latter shows that sufficiently far from the phase bound- 
ary t = the first order line is given by the relation 
v'^ = —6w^/t. For v' < Vc the effective action is min- 
imized by a configuration, where (j) = and ip points 
along one of the cube edges. This is the SSI* phase. 
For v' < Vc both the checkerboard and stripe domain 
order parameters become finite and ip points along one 
of the cube diagonals, as is the case for the SS2 phase. 
In particular one finds that at the phase transition from 
the superfluid to the SS2 phase ordering of the i/j-field 
induces ordering in the Ising field (j). In mean-field (mf) 
approximation one gets = {—w){ip)'^/t, due to the 
trilinear coupling to the ■0-field. 

The competition between the quartic and quadratic cu- 
bic symmetry breaking fields also leads to the existence 
of an intermediate supersolid phase (SS3). In the SSI* 
phase the quartic field w' < dominates and we have 
alignment along one of the cube edges in the {ipi,ip2)- 
planc. Right at the phase boundary between the SSI* 
and the SS3 phase the spin-wave spectrum becomes soft 
at the checkerboard ordering wave vector qq (i.e. t = 0). 
Hence within the SS3 phase we have a finite mean value 
for the Ising field, (0) 7^ 0, which implies that in addi- 
tion to the quartic cubic anisotropy {v' < 0) we have a 
quadratic term, gipi'ip2 with g — 'w{(j)). The latter favors 
spin alignment along the cube diagonals in the (^1, "02)- 
plane, which now starts competing with the quartic sym- 
metry breaking field v' < 0. By increasing the magnetic 
field and moving away from the SS1*-SS3 phase bound- 
ary the expectation value of the Ising field (j) and hence 
the magnitude of the quadratic symmetry field g starts to 
grow leading to a rotation of the spin alignment from the 
cube edges to the cube diagonals (compare also Fig. [l^ ). 
When the rotation is completed there is a phase tran- 



sition into the SS2 phase. With increasing U2/U1 the 
bicritical point on the tip of the SS3 phase meets with 
the phase boundary to the cantcd-ferromagnetic phase 
and it becomes a tetracritical point p^ . 

The above symmetry considerations allow us now to 
discuss the universality classes of the phase transitions 
between the various supersolid phases and the superfluid 
phase. Let us first discuss the phase transitions between 
the superfluid and the supersolid phases. The transition 
from the superfluid to the SSI* phase (coUinear super- 
solid = CSS) is in the universality class of the D — 3 
classical XY model ( CSS-SF universality class) Q. At 
the second order phase boundary between the superfluid 
and the SS2 phase (mixed supersolid = MSS) the checker- 
board and the stripe domain ordering wave vectors be- 
come soft simultaneously. As discussed above, this is not 
a consequence of both t and t' becoming zero at this phase 
boundary but due to the trilinear coupling between the 
checkerboard and stripe domain order parameters. Or- 
dering of the stripe domain fleld xp automatically induce 
ordering of the checkerboard field (see the above mean- 
field result). The critical exponent (3^ for the Ising field 
is (3^ — 1 already at the mean-field level. Whereas the 
critical behavior of the XF field seems to be unchanged 
at the transition, the trilinear coupling induces cusp-like 
singularities in the "slaved" checkerboard field (p. The 
critical behavior of such a model consitutes a new univer- 
slity class (MSS-SF universality class) and leave a more 
detailed analysis for a future investigation ||2^ . 

Next we consider the universality classes of the phase 
transitions between the various supersolid phases. By 
increasing the field h in the SSI* phase one may either 
encounter a first order phase boundary to the SS2 phase 
or a second order boundary to the SS3 phase. In the 
latter case the spin-wave spectrum becomes soft at the 
checkerboard ordering wave vector. Thus the effective ac- 
tion is given by an Ising-field which couples linearly to an 
XY field which is already ordered along one of the cube 
edges. After integrating out the ip field one is left with 
an Ising model with renormalized coefficients. Hence this 
transition belongs to the XSS-SF universality class 
In the intermediate supersolid SS3 phase the direction of 
the XY field tp rotates from the cube edge to the cube 
diagonal. The order parameter for the transition to the 
SS2 phase is the one-component field a — ipi ~ ip2- We 
thus conclude that this transition is again the XSS-SF 
universality class. 

At the multicritical point on the tip of the SS3 phase 
both of the Ising fields 4> and a = ipi — ip2 become criti- 
cal. The corresponding bicritical effective action for the 
charge ordering fields is given by 
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(3.8) 



Note that this is a rather pecuhar effective action for a 
bicritical point, since it consists of two Ising fields </> and 
(7 coupled by a three-point vertex. Since the upper crit- 
ical dimension of the coupling vertex wcjja'^ is Dc = 6 
as compared to Dc = 4 for the and vertices the 
former will dominate the critical behavior. In combina- 
tion with the coupling to the superfluid phase this should 
result in quite interesting critical behavior. The effec- 
tive action bears some resemblance with field theories 
for anisotropic Potts models [^. It would, however, go 
beyond the scope of the present paper to analyze the 
critical behavior of such a model. But, it certainly con- 
stitutes a quite interesting new quantum bicritical uni- 
versality class, which could be analyzed using standard 
renormalization group theory. 

With increasing the next-nearest neighbor interaction 
U2 the bicritical point moves towards the phase bound- 
ary to the superfluid phase. When it meets at about 
U2/U1 = 0.75 one gets a tetracritical point, where four 
second order lines meet. 

Finally, the phase transition from the supersolid SSI* 
phase to the stripe domain solid and from the SS2 phase 
to the 3/4 phase are both most likely mean-field like with 
logarithmic corrections in £) = 2 + 1 dimensions. 



C. Type III topology: Large next-nearest neighbor 
interaction (U2 > Ui) 

For this parameter range the zero-field classical ground 
state is an antiferromagnet with stripe symmetry. Due to 
the positive slope of the critical field in the SSI* phase 
(Eq. 2.39| ) the SS2 phase transition splits the discon- 
tinuous transition from the canted-ferromagnetic to 3/4 
phase into two continuous transitions. For U2 > 1.149 
this transition is established for the whole high field re- 
gion. 
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FIG. 13. Phase diagram for the hard-core Bose-Hubbard 
model with interaction parameters Ui = 1 and U2 ~ 1.5. The 
phase in between the SS2 and the SSI* is the SS3. 



(I) Supersolid phases and topology of the phase dia- 
gram: For values of the next-nearest neighbor interaction 
in the range of 1.149 < U2/U1 < 3 (see Fig. |3|) there are 
three different supersolid phases, SSI*, SS2, SS3. The 
SS3 phase is located in the thin sliver between the SS2 
and SSI* phase. Due to the presence of the intermedi- 
ate SS3 supersolid phases all phase transitions are con- 
tinuous. The transition from the SS2 phase directly to 
the canted-ferromagnetic phase without passing the SSI* 
phase has also quite interesting properties. Surprisingly, 
the transition line is given by 



U2-J 



U2 



(3.9) 



the same dependence as for the transition of canted- 
ferromagnet to SSI*. 

The topology for 1 < U2/U1 < 1.149 is similar to the 
just described one but with a slight change for the high 
field transition between the canted-ferromagnetic and the 
3/4 phase. When comparing the phase transiti on be - 
tween the canted-ferromagnetic and the 3/ 4 (E q. |2.26| ), 
the canted-ferromagne tic a nd the SS2 (Eq. |3.9| ) and the 
3/4 and the SS2 (Eq. |2.23D we encounter various region 
(as a function of J) where the SS2 phase occurs. In the 
region for U2/U1 < 1.124 we have for increasing J the fol- 
lowing picture (Fig. ^4|): There is a SS2 phase between 
the canted and the 3/4 (continuous on both transitions), 
then a SS2 phase with a discontinuous transition to the 
3/4 phase, a direct transition from canted-ferromagnetic 
to 3/4, a SS2 with discontinuous transition to the 3/4 
again and for high values of J the SS2 phase with con- 
tinuous transitions to either phases. 
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FIG. 14. Phase diagram for the hard-core Bose-Hubbard 
model with interaction parameters Ui — 1 and U2 = 1.1. 
The dashed transition line defines a discontinuous transi- 
tion from canted-ferromagnetic to 3/4, the solid lines for 
small J denotes two second order transition, namely from 
canted-ferromagnetic to SS2 and then to 3/4. The dotted lines 
denote a discontinuous transition from canted-ferromagnetic 
to SS2 and a continuous transition from SS2 to 3/4. The 
actual phase cannot be resolved on this scale. 
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The direct transition from canted- ferromagnetic to 3/4 
vanishes for U2/U1 > 1.124 at J = 0.415. For U2/U1 > 
1.149 only the SS2 phase with the continuous transitions 
keeps estabhshed. This hap pens for a criti cal va lue of 
J = 0.45. Equality of Eq. and Eq. define 

those multicritical points where the canted to 3/4 and 
the canted to SS2 phase lines meet. They are given by 
the two positive solutions of (besides J = 0) 

J3 + (2C/i - U2)J^ + Ui{Ui - 2C/2) J + t/2(C/| - Ul) = . 

For large values of next-nearest neighbor interaction, 
U2/U1 > 3.1, the phase diagram has a topology as de- 
picted in Fig. |l^ for the particular value of C/i = 1 and 
U2 = 5. The SS2 phase is now completely surrounded by 
the SSI* phase and there is no direct transition from the 
SS2 phase to the canted-ferromagnetic phase anymore. 
The tip of the SS2 lobe can be calculated analytically 
from the soft-mode of the SSI* phase (s. appendix 
It turns out that the value for Jtip has a remarkably sim- 
ple form, namely 



Jtip ~ Ul 



(3.10) 



In the extreme limit of vanishing nearest neighbor inter- 
action, Ul = 0, only the SSI* phase as a supersolid phase 
survives. The other supersohd phases and the 3/4 Mott 
insulator phase disappear (s. Fig. 

(11) Multicritical points and universality classes: The 
number of multicritical points in this parameter range is 
reduced. For U2/U1 < 3.1 there is a tetracritical point 
at the intersection point of the three supersolid phases 
(QC). This multicritical point disappears for U2/U1 ^3.1 
and there is a tricritical point at the tip of the SS2 phase 
instead. Furthermore there is a multicritical point at 
the tip of the Neel phase. The properties of all of these 
multicritical points have already been discussed in the 
preceding subsections. There are also no new universal- 
ity classes in the phase diagram which have not already 
been discussed. 
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FIG. 15. Phase diagram for the hard-core Bose-Hubbard 
model for U\ = \ and Ui = 5. The SS3 phase, which is lo- 
cated between the SS2 and the SSI* (lower fields) cannot be 
resolved at this scale. 
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FIG. 16. Phase diagram for the hard-core Bose-Hubbard 
model for pure next-nearest interaction, U2 = \ and U\ = 0. 
Note that there is only one supersolid phase (SSI*). 



IV. SUMMARY AND CONCLUSIONS 

The analysis of the hard-core Bose-Hubbard model 
with nearest and next-nearest neighbor interaction on a 
square lattice presented in this paper has demonstrated 
that the phase diagram can exhibit a rich variety of 
new quantum critical and multicritical phenomena. The 
topology of the phase diagram is determined by the in- 
terplay between the nearest and next-nearest charge in- 
teraction of the bosons and the cubic anisotropy. The 
principal effect of these terms is to stabilize several super- 
solid phases, which are characterized by the simultaneous 
presence of superfluidity and charge ordering, where the 
latter was described in terms of an Ising and a XY field. 

Topology of the phase diagram: We investigated the 
whole parameter range for nearest (Ui) and next-nearest 
neighbor (U2) interaction. By using ground state energy 
calculations and stability analysis (spin-wave theory) we 
have derived the phase diagrams in mean-field approxi- 
mation plus small fiuctuations. It turns out that there 
are three main types of phase diagram topologies. In pa- 
rameter regions {Ui > 2U2), where the Neel ground state 
is established at zero field (half-filling) , we encounter two 
supersolid phases, as already stated by Bruder et al. [D, a 
superfiuid phase, and three different types of Mott insu- 
lating phases. Whereas the spatial order of the supersolid 

551 phase can be described solely in terms of a checker- 
board pattern, the second supersolid phase (SS2) is a 
superposition of checkerboard and stripe domain density 
waves. In the region where a stripe ground state at half- 
fiUing is favored {U2 > Ui/2) three supersolid phases 
are found, SS2, SS3 and SSI*. The SSI* phase is the 
"supersolid partner" of the commensurate solid with a 
stripe domain structure and the spatial structure of its 
boson density can solely be described in terms of the 
stripe domain ordering wave vectors qi — 7r(l,0) and 
92 = 7r(0, 1). The phase transition from the SSI* to the 

552 phase is cither continuous or first order, depending 
on the particular choice of parameters. In terms of an 
effective action this can be restated as follows. The phase 
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transition from the SSI* to the SS2 phase is spht up 
into two phase transitions with an intervening supersoUd 
SS3 phase, if the Ising field describing the checkerboard 
structure orders before the first order transition occurs; 
otherwise the phase transition is first order. The tran- 
sition from the Mott insulator phase (3/4) is of second 
order into the SS2. For U2 < Ui the 3/4 phase joints 
partly (for higher fields) the canted-ferromagnetic phase; 
the transition is discontinuous. The occurrence of the su- 
persolid phases is strictly coupled to a finite value of U2 ■ 
Even in the limit of vanishing nearest neighbor interac- 
tion, Ui = 0, there is the SSI* supersolid phase present. 
It is also quite remarkable to note that in the parame- 
ter range 1.149 < U2/U1 < 3 the transition line from 
the canted ferromagnetic to the SS2 phase is the same as 
for the transition from the canted ferromagnetic phase to 
the SSI* phase over the whole range of hopping matrix 
elements J. 

Figs. |l^^ show the complex topology of the phase 
diagram for fixed hopping matrix elements but varying 
strength of the next-nearest neighbor interaction. 
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FIG. 17. Phase diagram for the hard-core Bose-Hubbard 
model with C/i = 1 and J = 0.4 as a function of the 
next-nearest neighbor interaction (72. 
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FIG. 18. The same as Fig. ^but zoomed. The SSI* phase 
is stable for slightly smaller values than (72 = 0.5, which marks 
the transition between the Neel and the stripe phase. 



In investigating the critical and multicritical behavior 
in the various parameter regions of the hard-core Bose- 
Hubbard model we have also found a rich variety of crit- 
ical and multicritical phenomena. 

Superfluid to supersolid transition: There are three dif- 
ferent universality classes for the phase transition from 
the superfluid to the supersolid phases. These are the 
XSS-SF, CSS-SF and MSS-SF universality class for the 
transition from the superfluid (SF) to the checkerboard 
supersolid (SSI), stripe domain supersoHd (SSI*) and su- 
persolid with an intermediate boson density profile (SS2). 
The CSS-SF transition is in the universality classes of 
the D = 2 + 1 = 3 dimensional XY model. The critical 
fluctuations do not lead to deviations from Bose-liquid 
behavior at the transition. The critical phenomena at 
the XSS-SF transition are more interesting since it can 
either be fluctuation driven first order or be governed 
by a strong coupling fixed point implying that the su- 
perfluid component displays "non Bose-liquid" behavior. 
For phase diagrams of type I topology the phase transi- 
tion from the superfluid to the mixed supersolid phase 
is always first order. Starting with U2 ~ Ui/2 (type II 
topology) one finds that within mean-field theory part of 
this phase boundary becomes second order and a tricrit- 
ical point appears. The critical behavior at the second 
order line constitutes a new MSS-SF universality class, 
which is interesting in several respects. First of all, the 
trilinear coupling between the checkerboard and stripe 
domain order parameter induces criticality even for a 
massive checkerboard field at the phase boundary where 
the mass of the stripe domain field becomes zero. Sec- 
ond, the critical exponent for the "slaved" checkerboard 
field 4> are anomalously large already at the mean-field 
level. For instance one finds 0nif oa (V')^, which implies 
/9™^ = 1. We suppose that a detailed analysis ||2^] of 
this new universality class will reveal interesting critical 
anomalies; in particular there will be cusp-like singulari- 
ties in the "slaved" checkerboard field. 

Supersolid to supersolid transition: The phase transi- 
tions between the various supersolid phases can either be 
first or second order. If they are second order, they all be- 
long to the XSS-SF universality class. In the parameter 
range where the checkerboard solid is the ground state at 
half-filling (type I topology) there is a first order transi- 
tion from the checkerboard supersolid (SSI) to the mixed 
supersolid phase (SS2) at large fields (fillings). Lowering 
the field the transition becomes continuous at the tip of 
the SS2 phase. The critical behavior at this point is de- 
scribed by a new quantum tricritical model different from 
the D = 3 classical Ising model. In the parameter range 
where the stripe domain solid is the ground state at half- 
filling there is no direct second order transition from the 
stripe domain supersolid (SSI*) to the mixed supersolid 
(SS2); it is always first order. At certain regions of pa- 
rameter space, however, there is an intermediate super- 
solid phase (SS3) and the first order line splits into two 
second order phase boundaries (both are in the XSS-SF 
universality class). At the tip of the intermediate SS3 
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phase there is a new bicritical quantum model. The ef- 
fective action contains two Ising fields which are coupled 
by an anisotropic trilinear vertex. The structure of the 
effective action shows some resemblance with anisotropic 
Potts models. 

Supersolid to commensurate solid transition: The 
phase transitions from the supersolid to the neighboring 
commensurate solid phases are all most likely mean-field 
like with logarithmic corrections in D = 2-1-1 dimensions. 

Multicritical points: For finite values of the hopping 
matrix element we have found up to seven different mul- 
ticritical points on a single phase diagram. Some of them 
are novel type of quantum multicritical points. In par- 
ticular, we found a new quantum tricritical point at the 
tip of the mixed supersolid phase. In the vicinity of these 
tricritical points there are critical end points with critical 
properties which are an interesting field for future inves- 
tigations. Especially the critical end point with a criti- 
cal "spectator phase" might engender novel "non Bose- 
liquid" behavior in the spectating phase. 
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APPENDIX A: SPIN- WAVE ANALYSIS OF THE 
3/4 PHASE 

In this section we present the main results of 
the spin- wave analysis of the 3/4 phase. Due 
to the four-sublattice structure we introduce four 
different Boson creation and annihilation operators 
(a)' , a; ) , (fel , 6; ) , (c| , Q ) , {dj ,di). By using the linearized 
Holstein-Primakoff transformation for the first three with 
spin up {I G Ci,C2 and £3) and the fourth spin down 
(m e £4), 



sr = 



Sf = 



Sf 

CZ 

QX 

qy 



-m ^ \ ai 
S + d\di 

2 

's 



{di+4) 
di - 4) , 



(Ala) 
(Alb) 

(Ale) 
(Aid) 

(Ale) 
(Alf) 



and inserting in the Hamiltonian (Eq. 1.2) we derive the 
following Hamiltonian: 



n 



NSh 



q 

+ Bq{aqbl + a^bq + Cqd-q + C^d'^-q) 
+ Cq{aqd-q + a^dl^q + bqC^ + b^Cq) 



with the coefficients 



8SU2 
^ h + 8S{U2~Ui) 
A^-^^ ^-h + 8S{Ui + U2) 

Bq = —AS J COS Qx 
Cq — —AS J COS Qy . 



(A2) 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 



This bilinear Hamiltonian can be diagonalized by stan- 
dard method, e.g. by introducing proper Greens func- 
tions. The excitation spectrum is then given by the zero 
of the following determinant 



V 





Bq 


Cq 




Bq 







Cq 


Cq 







LO Bq 





Cq 


Bq 


-UJ 1 



(A4) 



As a result, we get a polynomial of fourth order in the 
excitation energy 



7-a; 



tu! ■ 







(A5) 



with coefficients 



^2 



2.1(1 



(A6a) 

.1(1)' + 2A(i)A(2) - 2.1(1)^(4) _ .1(2)^(4) (A6b) 
2^(1)^(2)^(4) +A(i)'(yl(4)_A(2)) 

-(B2-C2)(A(2)+A(4)) 



(A6c) 



-A(1)'.1(2)A(4) +.1(1)(A(2) + A(4))(B2 + C2) 



-{B'q 



C: 



2\2 



(A6d) 



This polynomial defines the four branches of the exci- 
tation spectrum. Calculating the soft-mode at g = 
leads to the phase boundaries of the 3/4 phase. Due to 
Co = Bq the spectrum can be simplified: the first solu- 
tion is w = A^^\ independent of q. The other sol utions 
are given by a third order polynomial given by Eq. 2.23| . 
For J = the solution are given by ^(^^ , ^(2) and^4^ 
from which the boundaries can be readily seen. 



APPENDIX B: SPIN- WAVE ANALYSIS OF THE 
SSI AND SSI* PHASE 

These phases have two-sublattice spin structure which 
can be characterized by the wave vectors q which is go for 
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the SSI phase and qi and q2 for the SSI* phase. Spin- 
wave calculation is done via the linearized transformation 



in Fourier space (Eq. 1.6c). The biquadratic term in the 
Bose operators can be written in the form |24] 



For the SSI phase, q — qq, a soft-mode appears at qi and 
(72- The spectrum for this wave vector can be simplified 

(2) 

and vanishing of the second branch, Eq , leads to the 
condition 



Hi2) 



t t 

an a-n + al.a' 



'9 -9 



+ Cq (aq a-q-q + ajalg_^ 
+ Dq (ajflq+g + 



with coefficients 



2 Jo — Uq — Jq 



B„^S 



{2Jq + 2Jo - 2Uq - 2Uo 

Uq+q -Uq + Jq+q + Jq) SVO? 7 sin^ 5 
{2Uq - Jq+q - 2 Jq + Uq) Siu^ ,5 

(2C/o - 2 Jo -Jq + Uq) sin^ 7 
h sin 7 COS S 



Uq-Jq 



+ {Uq + Uq+q - Jq - Jq+qo) Slu^ 7 sin^ S 
+ {jq+q - Uq) slu^ S + {jq - Uq) siu^ 7 
Cq = S[jq — Uq) COS 7 siu (5 Slu 7 COS (5 
Dq = S\{Uq + Uo-Jq-Jo + Uq~ Jq) 



X COS 7 sin S sin 7 cos S 



-— cos 7 sin S. 
2 ' 



(B2c) 



Here the angles S and 7 are the difference and su m an gles 
of t he sp ins in respect to the magnetic field (Eq. 2.30 and 
Eq. 2.41 ). The excitation spectrum for this Hamiltonian 
then yields 



= ^{ni±n2)/2, 



(B3) 



with 
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A. 



91 



B. 
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2{Dq,-CqJ=0 



(B5) 



which defines with help of the angles 7 and S a cubic 
equation in the field h. The tip of the SS2 lobe then is 
given by that point where two solutions are degenerated, 
(Bl) i.e. where the discriminant vanishes. 

In the striped phase SSI*, q = 91,(72, there is a soft- 
mode at the checkerboard wave vector q — qo. In this 
case the spectrum becomes a quartic function in the field 
and the tip of the SS2 lobe can be found by studying the 
corresponding cubic resolvent. 
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